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Abstract 

We evaluate the potential of the ATLAS detector for discovering black holes 
produced at the LHC, as predicted in models with large extra dimensions where 
quantum gravity is at the TeV scale. We assume that black holes decay by Hawking 
evaporation to all Standard Model particles democratically. We comment on the 
possibility to estimate the Planck scale. 



1 Introduction 



Black holes will be produced at the Large Hadron Collider (LHC) if the fundamental 
Planck scale is of order a TeV p[] |2j. This scenario occurs in the model of large extra 
dimensions |Hj proposed by Arkani-Hamed, Dimopoulos and Dvali, where the Standard 
Model (SM) gauge and matter fields are confined to a 3-dimensional brane while gravity 
is free to propagate in extra dimensions of large size. This model was motivated by the 
necessity to solve the hierarchy problem. 

Using Gauss' law, the fundamental Planck scale of the (4+n)-dimension Mp is related 
to the 4-dimensional Planck scale Mpi (~ 10 19 GeV) by 

Mpj ~ M" +2 iT, 

where n is the number of extra dimensions and R is the size of the compactified dimen- 
sions. Assuming that the fundamental Planck scale Mp is the same as the electroweak 
scale (~1 TeV), the case n = 1 yields a very large R (~ 10 13 cm), ruled out by experi- 
ments. For n > 2, the size of R is less than ~ 10 -2 cm, which does not contradict results 
of gravitational experiments [3]. From astrophysical constraints j^j, the size of Mp is 
larger than O(TeV) and O(10 TeV) for n = 2 and 3, respectively. In particular, the 
case of n = 2 would be ruled out from the point of view of the solution of the hierarchy 
problem jE]. 

We consider here black holes of mass Mbh much larger than the fundamental Planck 
scale Mp since the phenomenology of black holes at Mbh ~ Mp is very complex and 
beyond the scope of the present study. 

For the collision of two partons at a center-of-mass energy = Mbh, by classical 
arguments £Q the total cross section is given by 

a (M BH ) ~ kR 2 s - 



Mp v M P K n 



l + n 



where R$ is the Schwarzschild radius: 



i r M BH ,8r(2p) 

R s = —= ( — ) 

b v^FMp 1 M P K n + 2 n 

For proton-proton collisions at the LHC, the differential cross section is given by 

dajpp -> BH + X) dL 
dM^ = dM^ a{ab ^ BH )I— Ih' 



'BH "iWBH 

dL 2M BH f 1 dx a , . , „ ,M 2 



BH 



dM BH s *—?Jm*i s x " sx « 

a.b BH»° 



a{ab^ BH)| 5=M | h = tiR, 



where a and b are partons in protons and fi(x) are the parton distribution functions (PDFs). 
An exponential suppression of the geometrical cross section has been proposed by Voloshin [7] . 
However since subsequent studies did not support this result [H], we do not consider the 
effect here. 
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Black holes decay through several phases [2j: balding, Hawking evaporation (spin-down 
and Schwarzschild) and Planck phases. Since the Hawking evaporation phase is expected 
to be the main phase of the decay, we naively consider only this process and assume that 
a black hole evaporates until its mass becomes zero. 

The radiation is characterized by the Hawking temperature T#, 

™ , , r w s , r ,M P n + 2 . i n + l n+1 
T H (M P ,n,M Bn ) = M P ( — -— o- ^ — F-" = • 

The heavier the black hole, the colder are its decay products. As it evaporates, its 
temperature increases but we will ignore the time evolution of black holes in this phase. 

We assume that the evaporation is described by black body radiation. The energy 
spectrum of the decay products obeys the Boltzmann distribution [T]: 

dN x 2 
dE ~ e x + c 

where x = E/Th and c is a constant, which depends on the quantum statistics of the decay 
products, i.e., c = — 1,+1,0 for bosons, fermions and Boltzmann statistics, respectively. 
In this study, we use c = for all particles. 

Certain conservation laws are obeyed in the decay of back holes. We assume that 
black holes decay "democratically", i.e., with roughly equal probability to all of the SM 
particles. 

In this paper, we present the discovery potential of the black holes at the LHC We 
describe a method of estimation of the fundamental Planck scale Mp from the discovery 
potential. 

The outline of this paper is as follows: the simulation conditions and description of 
signal and background samples are given in Section 121 Our original generator of the black 
holes is explained in detail in Sectional In Section EJ we describe the event selection and 
the reconstruction and finally, a conclusion is given in Section EJ 



2 Simulation 

The generator, developed for the purpose of this study, is described in Section El Initial 
state parton showers, hadronisation and decay are performed using PYTHIA 6.2 9J. All 
background samples are generated by PYTHIA 6.2. CTEQ5L is used for the parton 
distribution function. 

All samples are processed through the parameterized simulator, ATLFAST ^0], of 
the ATLAS detector. The energy resolutions and efficiencies for jet and particle recon- 
struction are corrected using the results of the full detector simulation. 

2.1 Signal 

We have generated signal samples for the various values of (Mp, n) listed in Tabled 
We have generated black holes whose mass Mbh is larger than Mp. Factorization and 
renormalization scales are set to a mass Mbh of the generated black hole. We call the 
lower limit of a black hole mass a threshold mass M B g = / B jj x Mp. We assume that 
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the energy spectrum of all the products from the black hole decay obeys the Boltzmann 
distribution. 



Table 1: Cross section of signal samples for each (Mp,n) parameter, f^, which is 
described in the text, is 1 for all cases. The unit of Mp is TeV. 



M P ,n a (pb) M P ,n a (pb) 



1,2 


9.45 • 10 3 


5,2 


0.662 


1,3 


8.26 • 10 3 


5,3 


0.603 


1,4 


8.06 • 10 3 


5,5 


0.625 


1,5 


8.24 • 10 3 


5,7 


0.699 


1,7 


9.05 ■ 10 3 






3,2 


25.5 


6,2 


0.125 


3,3 


22.9 


6,3 


0.114 


3,5 


23.4 


6,5 


0.119 


3,7 


26.0 


6,7 


0.133 


4,2 


3.74 


7,2 


0.0229 


4,3 


3.38 


7,3 


0.0210 


4,5 


3.49 


7,5 


0.0220 


4,7 


3.89 


7,7 


0.0247 



2.2 Background 

Our background samples are listed in Table El We use kinematical cuts at generation 
level to save CPU time and storage space of the data. The lower limit for the energy 
in the center-of-mass of initial partons, is set to 50 GeV for all the samples. The lower 
limit for the transverse momentum in the rest frame of initial partons pjf m , is also used 
as shown in Table El 

Table 2: Cross section of background samples for each mode. "Proc" is a process identi- 
fication defined in the PYTHIA. 



Process a (pb) p™ m Proc 



qq (q = 


=quark,lepton,gluon) 


1.29 • 10 4 


280.0 


11,12,13,28,53,68 




tt 


493 


10.0 


81,82 






0.468 


240.0 


25 




w ± z 


25.9 


10.0 


23 




zz 


10.6 


10.0 


22 




77 


229 


10.0 


18,114 




(V = w ± ,r,z) 


280 


10.0 


19,20 




W ± q 


73.4 


240.0 


16,31 




Zq,l*q 


31.5 


240.0 


15,30 




7? 


23.5 


240.0 


14,29,115 
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3 Generator 



A black hole generator was developed for this study, based on the assumptions and 
approximations discussed in Section QJ. 

We require the following conditions for the decay of black holes: 

• The number of fermions constrained by spin conservation 

• Four-Momentum conservation 

• Assumption of Boltzmann distribution for the energy spectrum of decay products 

• Color conservation 

• Assumption of democratic decay 

• Charge conservation 

• Option of conservation of the difference between the lepton number and the baryon 
number (B — L conservation) 

• Option of implementing Voloshin suppression [2j 

We can select whether to apply the B — L conservation and the Voloshin suppression, 
but they are switched off (no B — L conservation and no Voloshin suppression) in this 
study. 

For the production of the black holes, we use the Eq. ((TJ) of Section [T] Using the 
PDFs, the types of partons in the initial hard scattering are selected randomly. 

3.1 Spin Conservation 

We require that the number of fermions in the decay products of the black holes be even 
or odd according to the parity of the number of initial state fermions. Note that we do 
not consider gravitons in this study. 

3.2 Four-Momentum Conservation and Boltzmann Distribution 

We assume that the energy spectrum of the decay products follows a Boltzmann distri- 
bution. However, since energy and momentum must be conserved, all particles cannot 
be sampled randomly from that distribution. We do not allow the case N = 2 since 
it is fully constrained. For N > 3, particle energies are successively sampled from the 
Boltzmann distribution until the summed energy exceeds the mass of the black hole. The 
last particle is then given the energy required by the energy constraint. The first N — 2 
particles are assigned random directions and the last two are emitted in a direction such 
as to conserve overall momentum. 

1 Assumptions and approximations of our generator are similar with a generator, TRUENOIR, de- 
scribed in Another generator has recently been published \12\ . 
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3.3 Color Conservation 



Color is assumed conserved in the decay of black holes. Color is assigned to the initial 
partons and color connection is implemented. As an example, if the initial partons are 
two gluons, we select randomly either a qq pair or a gluon from the final state partons, 
if present. We then connect the colors as shown in Figure H If no qq pair or gluon is 
present, the color connection must be applied between the two initial state gluons, if 
possible. Another example of color connection when the initial partons are a quark and 
a gluon is shown in Figure ^ 

For the remaining partons, color connection is applied between pairs with closest 
opening angle. If there is an odd gluon remaining, it is connected to the qq pair or a gg 
pairs or to the initial partons. 




Figure 1: Color connection for gg (left) and gq (right). 



3.4 Democratic Decay 

We assume a democratic decay, constrained by the conservation laws. The number of 
degrees of freedom for each particle takes into account charge, spin and color. We do not 
consider gravitons in this study. The mass of the SM Higgs is set at 120 GeV. Events are 
produced where the types of particles are chosen randomly with the probabilities listed 
in Table El and are accepted only if conservation laws can be applied. 
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Table 3: Degrees of freedom and assigned probability in the generator for each particle. 



Particle Degrees of freedom Assigned probability 



g (gluon) 


8 


0.0690 


W 


(3 


0.0517 


Z 


3 


0.0259 


7 


2 


0.0172 


lepton (e, /i, r) 


4 


0.0345 


neutrino (v e , z/ M , u T ) 


4 


0.0345 


quark (u, d, c, s, t, b) 


12 


0.1034 


Higgs 


1 


0.0086 


Graviton 


5 


0.0000 



3.5 Performance 

Figure |2 shows the differential cross sections, calculated from Eq. (JjQ). Figures H3 show 
the shape of the mass distributions of generated black holes. The distributions shown 
start at Mbh > Mp although, as mentioned above, the validity of the model applies in 
the region M B h 3> M P . 



> 




M BH (TeV) 



Figure 2: Differential cross section as a function of a black hole mass for each (Mp,n) 
parameter. 

Figures |U show various distributions of the generated black holes: p z , charge, lepton 
number (L), baryon number (B), B — L, and the multiplicity of decay products. The 
absolute values of L are always even because of the spin conservation as shown at Fig- 
ures IU(c) and (g). When spin conservation is turned off, they take on both even and odd 
values, as shown at Figures |U (h). Figures |U (f),(i) and (j) indicate that the multiplicity 
of decay products depends on M B h not M P . 
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(a) (M P , n)=(l,3) 



(b) (M P , n)=(l,7) 




M BH (TeV) M BH (TeV) 

(c) (Mp, n)=(3,3) (d) (M P , n)=(5,3) 




Figure 3: Mass distributions of the generated black holes. The histogram shows generated 
black holes and the solid line shows an expected shape from Eq. (|T|) . 

The distribution of particle types in the decay of black holes are shown in Figures El 
As expected, there are more quarks than anti-quarks (Figures El (a) and (b)) since the 
LHC is a proton-proton collider. We have also generated signal samples with (Mp, n, 
/bh) = (1'3,9) in order to analyze events with high multiplicity of decay products. As 
the multiplicity of decay products becomes larger, the proportion of the different types 
of products becomes closer to the values of Table El even when many conservations are 
imposed, as can be seen shown in Figures El (c) and (d). When there is no charge and 
spin conservation, the ratios are almost the same as the values of Table El as shown at 
Figures El (e) and (f). 
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(a) p z for (1,3,1) 
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(j) Multiplicity for (7,3,1) 
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Figure 4: Distributions of generated black hole properties: p z , charge, lepton number (L), 
baryon number (B), B — L, and multiplicity of decay products, (a)-(f) for (Mp, n, 
/* h H ) = (l,3,l), (g)-(i) for (M P , n, /^ H ) = (l,3,9) and (j) for (M P , n, ^ H ) = (7,3,l). 
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(a) PDG code for (1,3,1) 



(b) |PDG code| for (1,3,1) 
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(c) PDG code for (1,3,9) 
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(d) |PDG code| for (1,3,9) 
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Figure 5: Distributions of PDG code of decay products of black holes. The left shows 
PDG code and the right shows absolute value of PDG code. The points of * show the 
values at Table El 
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Figures El show various distributions of the decay products: px, p z , energy, 77 and 4>. 
Distributions of p z , 77 and of the products are symmetric as expected. 



(a) p T for (1,3,1) 
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Figure 6: Distributions of various properties of decay products of black holes: pr, Vz-, 
energy, rj and (j) 



4 Analysis 

Although the expressions for black hole production given in Section ^ are not valid at 
Mbh ~ Mp, owe BH generator applies them in all the regions where Mbh > Mp in this 
study. 
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4.1 Selection Criteria and Reconstruction 

The following cuts are applied in order to select and reconstruct black holes (BHs): 

• For a precise reconstruction of the BH, it is necessary to remove particles produced 
in the stage of initial state radiations (ISRs). Because the particles from ISRs tend 
to have small px and large \r]\, as shown in Figures. [7|and|SJ the following require- 
ments are applied to particles in each event (ISR-cut) : 

• pr > 30 GeV for /i, e 

• pt > 50 GeV for 7, jet 




20 40 60 80 100120140160180200 
Pt of p.,e(GeV) 
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CD 
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E 

o 
c 



10 



10 



-2 



solid - from BH 

(Mp=1TeV,n=3) 
dashed - from 

initial radiation 



20 40 60 80 1 00120140160180200 
Pt of yJets(GeV) 



Figure 7: Pt distributions : comparison between particles from BH and initial state 
radiation (Mp = 1 TeV, n=3). 
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Figure 8: rj distributions : comparison between particles from BH and initial state radi- 
ation (M P = 1 TeV, n=3). 

• Of all the particles passing the ISR-cut, more than three are required to have energy 
larger than 300 GeV (E > 300 GeV), and moreover, at least one of them has to be 
either an electron or a photon. This latter requirement is to suppress backgrounds. 
Figure El shows the distributions of the number of energetic particles in each event, 
where event is one of the largest contribution in the backgrounds as shown in 
Table |U It is seen that this multiplicity cut on energetic particles is very effective. 

• We require the event shape variable R 2 to be less than 0.8 (R 2 < 0.8). R 2 represents 
an event topology and it is defined by Fox- Wolfram moments as follows: 

R 2 = H 2 /H 

where Hi is the i-th Fox- Wolfram moments, j and k are ID numbers of tracks, p* 
is the momentum of the track j in the rest frame of BH, <pjk is the opening angle 
between the tracks j and k , and Pi(x) is the Legendre polynomial. E is obtained 
by summing up the energies of particles passing the ISR-cut, calculated in the rest 
frame of BH. R 2 ranges from to 1 (0 < R 2 < 1) as can be seen from Figure fTUl 
which shows the R 2 distribution for signal and qq events with or without selection 
criteria. Since lower value of R 2 indicates more spherical event, we remove qq events 
here. 
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Figure 9: Shape of the distribution of the number of energetic particles with E > 
300 GeV : comparison between signal and 75 background (M P = 1 TeV, n—3). 




Figure 10: R2 Distributions : comparison between signal event and qq event (Mp = 1 TeV, 
n=3). Left: without any selection criteria and right: with selection criteria 

• Et < 100 GeV. For a precise calculation of BH mass, events with high missing 
energy are rejected here. 
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The black hole mass is then reconstructed from the 4-momenta of the remaining 
muons, electrons, gammas and jets as follows: 

Per = Pi 

i=/x,e,7,jet 
M BH = SJP^H, 

where Pi is a reconstructed four-momentum of each particle. 

The quality of the reconstructed mass of the BH is discussed in Appendix [X] 

4.2 Discovery Potential 

We have evaluated the discovery potential for the following cases of (Mp ,n): 

M P = l,3,4,5,6,7TeV 
n = 2,3,5,7 

by using the selection criteria described in Section 14.11 Because Mp is an unknown 
parameter, we cannot set the proper cut value for the lower limit of Mbh- Therefore we 
consider various values as follows and evaluate S/ \f~B in each case: 

M BH >M^ n = l,2,3,4,5TeV. 

The Mbh distributions are shown in Figure ITT1 for different values of Mp with n = 3. 
These are obtained after all the selection criteria have been applied. The histograms 
represent the sum of signals and backgrounds with the cross-hatched part showing only 
backgrounds. Note that Mgjj 11 is set to 1 TeV for this figure. 
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(a) M P =1 TeV 



(b) M P =3 TeV 





Figure 11: Mbh distributions in the case of n=3 and M^=l TeV. (solid line : signal 
plus background, cross hatched : background only) 
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Table 4: The number of remaining BG events at J Cdt=10 fb 1 . e M ™ 1 is the value of 
efficiency for M™ n = 1 TeV. 





qq 


tt 






a (pb) 


1.29 • 10 4 


493 




0.468 


M mm 
£ JU BH 


=i 


1.67- 10~ 5 


1.39-10 




8.5 • 10~ 5 




1.0 


2.15 ■ 10 3 


68.5 




0.40 




2.0 


1.18 ■ 10 3 


13 




0.1 




3.0 


4.3 • 10 2 


2 




0.02 


(TeV) 


4.0 


77 


0.5 









5.0 


3- 10 














W ± Z 


ZZ 


77 


7 y 


a (pb) 


25.9 


10.6 


229 


280 


ll,fmin 
£ JlJ BH 


=i 


5-10~ 7 1.5 • lO" 6 





7- 10" 7 




1.0 


0.1 


0.084 





2 




2.0 


0.05 


0.04 





1 




3.0 














(TeV) 


4.0 
















5.0 

















W ± q 


Zq,l*q 


1Q 


total 


(pb) 


73.4 


31.5 


23.5 






4.7- 10~ 5 


1.40 • 10~ 4 


1.484- 10~ 3 




1.0 


35 


44.2 


349.3 


2.66 ■ 10 3 


M™ 11 2.0 


18 


21 


195 


1.43 • 10 3 


3.0 


4 


3.5 


49.2 


4.9 ■ 10 2 


(TeV) 4.0 


0.7 


0.3 


7.8 


86 


5.0 








2.8 


3-10 



Tables EE1 show the number of events passing the selection criteria for background (BG) 
processes and signal events with various (Mp ,n). The number is normalized for J Cdt=10 fb 
The values of efficiency are listed only for M^" = 1 TeV. As can be seen in Table EJ 
7g process has the highest efficiency, but the background is dominated by the qq process 
owing to its large cross section. The total number of BG events is found to be ~ 10 3 - 10 4 
in the case of Mgf 1 = 1-2 TeV and ~ 10 - 10 2 in the case of = 4-5 TeV. 

According to Tables El and |H1 the efficiency for signal events is ~ 1% in all cases of 
{Mp ,n). It is relatively lower when Mp = 1 TeV or 7 TeV because signal events are 
more likely to be rejected by the M^ cut in the case of smaller Mp, or they have a large 
I$t f° r larger Mp. 

We calculate S/ \/~B from the data of Tables IU- El and estimate J Cdi BCO verydt — the 
integrated luminosity with which discovery is achieved. Here the condition of discovery 



16 



Table 5: The number of remaining signal events at j Cdt=10 fb 1 . £ m bh 1 is the value 
of efficiency for M™ n = 1 TeV. 



(M P ,n) 


(1,2) 


(3,2) 


(4,2) 


a (pb) 


9450 


25.51 


3.738 




9.63 x 10~ 3 


2.747 x 10~ 2 


2.733 x 10- 2 


(TeV) 


1.0 


9.10 x 10 5 


7008 


1022 


2.0 


6.50 x 10 5 


7002 


1021 


3.0 


2.52 x 10 b 


6594 


1016 


4.0 


7.5 x 10 4 


2987 


943.5 


5.0 


2.3 x 10 4 


1.06 x 10 3 


404.5 




(M P , n) 


(5,2) 


(6,2) 


(7,2) 




0.6622 


0.1249 


2.292 x 10~ 2 


£ M§™=1 


2.601 x 10~ 2 


2.287 x 10~ 2 


2.022 x 10" 2 


(TeV) 


1.0 


172.2 


28.56 


4.634 


2.0 


172.2 


28.54 


4.634 


3.0 


171.7 


28.53 


4.632 


4.0 


170.1 


28.46 


4.616 


5.0 


154.2 


28.15 


4.593 




(M P , n) 


(1,3) 


(3,3) 


(4,3) 


a (pb) 


8256 


22.89 


3.381 




8.46 x 10~ 3 


2.725 x lO" 2 


2.605 x 10^ 2 


(TeV) 


1.0 


6.98 x 10 5 


6238 


880.8 


2.0 


4.82 x 10 5 


6224 


879.4 


3.0 


2.01 x 10 5 


5878 


874.0 


4.0 


5.8 x 10 4 


2745 


810.8 


5.0 


1.2 x 10 4 


929 


354.3 




(M P ,n) 


(5,3) 


(6,3) 


(7,3) 


o (pb) 


0.6027 


0.1142 


2.105 x 10^ 2 




2.574 x 10~ 2 


2.347 x 10~ 2 


2.071 x 10" 2 


(TeV) 


1.0 


155.1 


26.80 


4.359 


2.0 


155.1 


26.80 


4.357 


3.0 


154.4 


26.76 


4.351 


4.0 


152.7 


26.59 


4.341 


5.0 


140.4 


26.23 


4.317 
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Table 6: The number of remaining signal events at j Cdt=10 fb 1 . £ m bh 1 is the value 
of efficiency for M™ n = 1 TeV. 



(M P ,n) 


(1,5) 


(3,5) 


(4,5) 


a (pb) 


8244 


23.42 


3.487 




7.24 x 10~ 3 


2.563 x 10~ 2 


2.779 x lO" 2 


(TeV) 


1.0 


5.97 x 10 5 


6003 


969.0 


2.0 


4.24 x 10 5 


5993 


969.0 


3.0 


1.72 x 10 b 


5600 


964.5 


4.0 


5.0 x 10 4 


2497 


892.0 


5.0 


1.4 x 10 4 


775 


377.3 




(M P , n) 


(5,5) 


(6,5) 


(7,5) 




0.6252 


0.1191 


2.203 x 10~ 2 


£ M§™=1 


2.501 x 10- 2 


2.244 x lO- 2 


2.139 x 10" 2 


(TeV) 


1.0 


156.4 


26.73 


4.712 


2.0 


156.3 


26.73 


4.712 


3.0 


155.7 


26.69 


4.706 


4.0 


154.9 


26.61 


4.688 


5.0 


139.2 


26.33 


4.662 




(M P , n) 


(1,7) 


(3,7) 


(4,7) 


a (pb) 


9053 


26.02 


3.887 




6.54 x 10~ 3 


2.577 x 10- 2 


2.630 x 10^ 2 


(TeV) 


1.0 


5.92 x 10 5 


6705 


1022 


2.0 


4.22 x 10 5 


6698 


1020 


3.0 


1.57 x 10 5 


6310 


1017 


4.0 


4.5 x 10 4 


2789 


943.0 


5.0 


5 x 10 3 


892 


382 




(M P ,n) 


(5,7) 


(6,7) 


(7,7) 


o (pb) 


0.6991 


0.1335 


2.474 x 10^ 2 




2.577 x lO- 2 


2.319 x 10~ 2 


2.090 x 10" 2 


(TeV) 


1.0 


180.2 


30.96 


5.171 


2.0 


180.2 


30.93 


5.171 


3.0 


179.7 


30.89 


5.168 


4.0 


178.3 


30.72 


5.151 


5.0 


164.7 


30.48 


5.124 
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Table 7: The value of S/VB at / Cdt=10 fb _1 . 



(M P , 


n) 


(1,2) 


(3,2) 


(4,2) 


(T (pb) 


9450 


25.51 


3.738 




1.0 


1.76 x 10 4 


136 


19.8 


M™ n 


2.0 


1.72 x 10 4 


185 


27.0 




3.0 


1.1 x 10 4 


3.0 x 10 2 


46 


(TeV) 


4.0 


8.1 x 10 3 


3.2 x 10 2 


1.0 x 10 2 




5.0 


4 x 10 3 


2 x 10 2 


7 x 10 



(M P , 


n) 


(5,2) 


(6,2) 


(7,2) 


a (pb) 


0.6622 


0.1249 


2.292 x 10~ 2 




1.0 


3.34 


0.554 


0.0899 




2.0 


4.55 


0.755 


0.123 




3.0 


7.8 


1.3 


0.21 


(TeV) 


4.0 


18 


3.1 


0.50 




5.0 


3 x 10 


5 


0.8 



(M P , 


n) 


(1,3) 


(3,3) 


(4,3) 


a (pb) 


8256 


22.89 


3.381 




1.0 


1.35 x 10 4 


121 


17.1 




2.0 


1.27 x 10 4 


165 


23.3 




3.0 


9.1 x 10 3 


2.7 x 10 2 


39 


(TeV) 


4.0 


6.3 x 10 3 


3.0 x 10 2 


87 




5.0 


2 x 10 3 


2 x 10 2 


6 x 10 



(Mp, 


n) 


(5,3) 


(6,3) 


(7,3) 


a (pb) 


0.6027 


0.1142 


2.105 x 10~ 2 




1.0 


3.01 


0.520 


0.0845 




2.0 


4.10 


0.709 


0.115 




3.0 


7.0 


1.2 


0.20 


(TeV) 


4.0 


16 


2.9 


0.47 




5.0 


3 x 10 


5 


0.8 
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Table 8: The value of S/VB at / Cdt=VQ fb -1 . 



(M P , 


n) 


(1,5) 


(3,5) 


(4,5) 


a (pb) 


8244 


23.42 


3.487 




1.0 


1.16 x 10 4 


116 


18.8 


M- n 


2.0 


1.12 x 10 4 


158 


25.6 




3.0 


7.8 x 10 3 


2.5 x 10 2 


44 


(TeV) 


4.0 


5.4 x 10 3 


2.7 x 10 2 


96 




5.0 


3 x 10 3 


1 x 10 2 


7 x 10 



(Mp, 


n) 


(5,5) 


(6,5) 


(7,5) 


a (pb) 


0.6252 


0.1191 


2.203 x 10~ 2 




1.0 


3.03 


0.518 


0.0914 




2.0 


4.13 


0.707 


0.125 




3.0 


7.0 


1.2 


0.21 


(TeV) 


4.0 


17 


2.9 


0.51 




5.0 


3 x 10 


5 


0.9 



(Mp, 


n) 


(1,7) 


(3,7) 


(4,7) 


<r (pb) 


9053 


26.02 


3.887 




1.0 


1.15 x 10 4 


130 


19.8 


M™ n 


2.0 


1.11 x 10 4 


177 


27.0 




3.0 


7.1 x 10 3 


2.9 x 10 2 


46 


(TeV) 


4.0 


4.9 x 10 3 


3.0 x 10 2 


1.0 x 10 2 




5.0 


9 x 10 2 


2 x 10 2 


7 x 10 



(Mp, 


n) 


(5,7) 


(6,7) 


(7,7) 


a (pb) 


0.6991 


0.1335 


2.474 x 10~ 2 




1.0 


3.49 


0.600 


0.100 


M- n 


2.0 


4.77 


0.818 


0.137 




3.0 


8.1 


1.4 


0.23 


(TeV) 


4.0 


19 


3.3 


0.56 




5.0 


3 x 10 


6 


0.9 
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is set as: 



S/VB > 5.0 and S > 10 



which is a conventional condition, as used in the analysis of Higgs events with ATLAS. 
The results for S/ VB and J C^i acoveIy dt are shown in Table 0|H] and El respectively. We 
require that M^ n be larger than Mp to calculate integrated luminosities for BH discovery. 
In Table the shaded values indicate the most favorable cut in M^ n in each case of 
(Mp ,n). From these tables, we see that the discovery can be accomplished within 
/ £discover y G^ < 1 fb" 1 in all cases of n if Mp is less than ~ 5 TeV. 

Figure H21 gives a contour plot for J Cd iscovcry dt in (M P ,n) plane. We find that the 
discovery potential hardly depends on n but has a strong dependence on M P . This is due 
to the fact that the cross section is a strong function of Mp but not of n, as shown in 
Tables El and El From this relation between Mp and J Cn scover ydt, Mp can be determined 
by how early the discovery is accomplished. 



> 

0) 



5.5 



1 fb" 



4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 



00pb~ 



0pb~ 



123456789 



Figure 12: Contours of J Cdincoverydt in (Mp, n) plane. 

In a quantitative point of view, it is found that the excess of events is detected in ~ 
1 month at low luminosity (J £dt=l fb -1 ) if Mp < 5 TeV, and discovery within only one 
day (/ Cdt=100 pb" 1 ) can be expected if M P < 4 TeV. 

As was previously mentioned, the BH model we assumed here is valid only when 
Mbh 3> Mp. As Mbh approaches Mp, the theory of BH production becomes very 
complex. If we consider events with reconstructed Mbh > 5 TeV in the case of Mp = 
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Table 9: Integrated luminosity (fb _1 ) required for discovery, for various values of (Mp, n), 
as a function of M^ n (M^ n > Mp). The shaded entry gives the most favorable value in 
each case of (M P ,n). 



(M P , 


n) 


(1,2) 


(3,2) 


(4,2) 


(5,2) 


<r(pb) 


9450 


25.51 


3.738 


0.6622 




1.0 


1.10 x 10" 4 










2.0 


1.54 x 10" 4 










3.0 


3.97 x 10~ 4 


1.517 x 10- 


-2 






(TeV) 


4.0 


1.3 x 10" 3 


3.348 x 10 


-2 


0.1060 






5.0 


4.3 x lO- 3 


9.43 x 10" 


2 


0.2472 


0.6485 




(M P , 


n) 


(1,3) 


(3,3) 


(4,3) 


(5,3) 


<r(pb) 


8256 


22.89 


3.381 


0.6027 




1.0 


1.43 x 10" 4 








M- n 


2.0 


2.07 x 10~ 4 










3.0 


4.98 x 10~ 4 


1.701 x 10- 


-2 






(TeV) 


4.0 


1.7 x 10" 3 


3.643 x 10- 


-2 


0.1233 






5.0 


8.3 x lO" 3 


0.108 


0.2822 


0.7123 



(M P , 


n) 


(1,5) 


(3,5) 


(4,5) 


(5,5) 


a (pb) 


8244 


23.42 


3.487 


0.6252 




1.0 


1.68 x 10" 4 








M^t 


2.0 


2.36 x 10~ 4 










3.0 


5.81 x 10~ 4 


1.786 x 10" 


-2 






(TeV) 


4.0 


2.0 x 10~ 3 


4.005 x 10' 


-2 


0.1121 






5.0 


7.1 x 10" 3 


0.129 


0.2650 


0.7184 




(Mp, 


n) 


(1,7) 


(3,7) 


(4,7) 


(5,7) 


a (pb) 


9053 


26.02 


3.887 


0.6991 




1.0 


1.69 x 10" 4 








M- n 


2.0 


2.37 x 10~ 4 










3.0 


6.37 x 10~ 4 


1.585 x 10" 


-2 






(TeV) 


4.0 


2.2 x 10~ 3 


3.586 x 10" 


-2 


0.1060 






5.0 


2 x 10" 2 


0.112 


0.262 


0.6072 
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1 TeV, we find that we need ~5 pb _1 instead of ~0.1 pb" 1 to discover black holes. 
However the events include black holes whose generated mass was less then 5 TeV. When 
black holes with mass greater than 5 TeV are generated, we find that we need ~10 pb" 1 . 

Considering that the model of BH formation and decay is valid only for Mbh 3> Mp, 
we show, in Figure EE an evaluation of the discovery potential of BHs when the cut 
Mg 1 ^ 11 > Mp + 1 TeV is applied. Although more integrated luminosity is required, it is 
clear that an excess of events will still be easily observed within a few days of running, 
for Mp values up to a few TeV. 



> 4 

0) 

^3.5 
3 
2.5 
2 



1.5 
1 



123456789 

n 



Figure 13: Contours of J £discover y <^ in (Mp, n) plane in case of M^ 11 > M P + 1 TeV. 

5 Conclusion 

We have studied the potential to observe black holes with the ATLAS detector at the 
LHC. We developed a generator of black holes with simple assumptions for the production 
and decay process, taking the expressions in Section Q as valid in the full region of the 
black hole mass Mbh > Mp. We find an excess of events from the Standard Model 
backgrounds for the integrated luminosity of 1 fb" 1 (~ 1 month at low luminosity) if 
Mp < 5 TeV and with the integrated luminosity of 100 pb _1 (~ 1 day) if Mp < 4 TeV. 
If we assume validity when Mbh > Mp + 1 TeV, the required luminosities are slightly 
larger. 
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We also proposed a simple method to estimate the Planck scale Mp. It will be nec- 
essary to find a way to estimate the number of large extra dimensions n and to identify 
the excess as the black holes. The development of more realistic black hole generators is 
also important for more detail studies. 
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A Reconstructed Mass of Black Hole 

The resolution in the reconstructed mass of the BH depends on various factors. We 
investigate reasons of the overestimation and the underestimation on the reconstructed 
mass Mbh as follows. 

Figures d show distributions of M BH - Mg£ c for BH events with M^ c > 1.0,6.0, 
and 9.0 TeV in case of (Mp,n) = (1,3), respectively. 

(a) (M P , n, A) = (l,3,l) (b) (M P , n, / B h H ) = (l,3,6) (c) (M P , n, /* h H ) = (l,3,9) 




Figure 14: Distributions of the difference between a reconstructed and a generated mass 
of a black hole Mbh — ^BtT- O ur selection criteria is applied to reconstruct a black hole. 
E miss in the figure means J^t- 

There are contaminations from particles arising from the initial state radiation and 
from the spectator proton fragments of the hard scattering. We call these particles 
"ISRs/SPECs" . FigureslTHlshow Mbh - Mg^ distributions of the events with no neutrino. 
We compare the difference in mass when all particles are used and when ISRs/SPECs 
are excluded. The generator information is used to determine whether a particle comes 
from ISRs/SPECs. We can see that the overestimation of a black hole mass is due to the 
contaminations of ISRs/SPECs. 
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a) (M P , n, /^ h H )=(l,3,l) (b) (M P , n, ^)=(1,3,6) (c) 



(M P , n, /&) 



;i,3,9) 



^ ......... -T 


Mean 306.1 

RMS 537.5 




-4000 -2000 2000 4000 

M BH -M BH fue (no neutrino) all 


! J 


Mean -26.77 
RMS 444.9 
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0.175 
0.15 
0.125 
0.1 
0.075 
0.05 
0.025 




0.175 
0.15 

0.125 
0.1 

0.075 
0.05 

0.025 



: J 


Mean 211.0 
. RMS 452.3 




-4000 -2000 2000 4000 

M BH -M BH tme (no neutrino} all 




RMS 509.3 






0.175 
0.15 

0.125 
0.1 

0.075 
0.05 

0.025 



-4000 -2000 2000 4000 

M BH -M BH ttue (no neutrino) all 
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M BH" M BH fme ( no neutrino) no ISRs/SPECs 



-4000 -2000 2000 4000 

M BH -M BH irue (no neutrino) no ISRs/SPECs 



-4000 -2000 2000 4000 

M BH" M BH lme ( no neutrino) no ISRs/SPECs 



Figure 15: Mbh - Mbk Distributions. We use events with no neutrino. The upper figure 
is obtained from all reconstructed particles and jets and the lower is from them except 
for ISRs/SPECs. Details are given in the text. 



Mgjf and a measured missing 



Figure shows the correlation between a Mbh 
energy We can see clearly that we underestimate a black hole mass if the missin 
energy is important. 



(a) (M P , n, /&)=(1,3,1) (b) (M P , n, /^ H )=(l,3,6) (c 
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Figure 16: Measured missing energy^ versus Mbh — ^biT- mass °f a black hole is 
obtained from reconstructed particles and jets except for ISRs/SPECs. 

We perform various corrections: (1) we subtract the momentum of particles which are 
selected by our criteria but come from ISRs/SPECs, from the reconstructed mass Mbh, 
(2) add the momentum of neutrinos, (3) add the momentum of particles which do not 
come from ISRs/SPECs but are removed by our selection criteria, as shown at Fig- 
ures El We can see that the reconstructed mass is improved by these corrections as we 
expected. From these results we understand that the overestimation of Mbh is caused by 
the contamination of particles from the initial state radiation and the spectators while 
the underestimation is due to the missing energy of neutrinos. 
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a) (M P , n, ^ h H )=(l,3,l) (b) (M P , n, /* h H )=(l,3,6) (c) (M P , n, /* h H )=(l,3,9) 
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0.05 



10 20( 



-2000 2000 
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Mean 25.44 
RMS 5043 



Mean -242.2 
HMS 5751 



Mean -425.S 
HMS MS 



0.05 




RMS 777.0 



)0 20( 



Mean -71 0.5 

RMS 7674 



M^-M^iE^O.ITeVl-ISRs'EPECs 



-2000 2000 

M SH -iM SH "(E rai!S <0.1 TeV) -ISRs/SPECs 





Mean 10.00 




RMS 495.0 


, I , , , I , tJ 5 


V™. , 1 i , , 1 i 



-4000 -2000 2000 < 

m bh" m bh i,ljb (E ni3S <0.1 TeV) -ISRs/SPECa+Neutrinos 



1000 -2000 2000 4000 

M BH" M BH tpje tW** 1 TeV ' - 1 SRs/SPECs+NeuhiFi os+OutOf Se lecbon 



Mean -258,3 
HMS 503.0 



-4000 -2000 2000 < 

Wti' m fi*s<M TeV ) -ISRs/SPECs+Neulrinos 



0.2 ; 
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Figure 17: Mbh — A^bh 6 Distributions. The top figures are the same with Figures IT41 
The 2nd, 3rd and bottom figures show M B h — distributions corrected step-by-step. 

Details are given in the text. 
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